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^ : Abstract 
O 

We present an on-shell renormalization scheme for the neutralino and chargino mass matrices 
for two i?-parity violating models, namely bilinear i?-parity violation and the /xz/SSM with 
. one right-handed neutrino superfield. We discuss in both models how to obtain neutrino 
masses and mixing angles correctly as well as differences to the existing DR calculations. 
^ ' Both models predict correlations between neutrino mixing angles and ratios of i?-parity 
violating decays such as xl ~^ I W T which can be tested for example at the LHC. We 
point out that there are cases where one has to include the complete NLO corrections for 
the decays xW ~^ l^W T to predict reliably the correlations between neutrino- and collider- 
physics. Moreover, we discuss briefly the implementation of these decays in the public 
CN ■ program CNNDecays. 
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1. Introduction 

The Minimal Supersymmetric extension of the Standard Model (MSSM) [l| is based on 
conserved i?-parity, where .R-parity (R p ) j2] is defined as R p = (— i} 3B + L+2s an d originally 
guaranteed the stability of the proton in supersymmetric models It implies on the 
one hand a stable lightest supersymmetric particle (LSP) and predicts on the other hand 
zero neutrino masses as in the Standard Model (SM). As neutrino oscillation experiments 
have shown at least two neutrinos have non-zero mass 0] requiring an extension to explain 
neutrino data0. 

In the SM exist various ways to explain neutrino data, e.g. postulating very heavy par- 



ticles as in the seesaw mechanism |ll|, |12| resulting in a unique dimension-5 operator [13 



or via loop-effects (12|, |14|, [15|. Although all mentioned solutions can be supersymmetrized, 



there exists an entirely supersymmetric possibility to generate Majorana neutrino masses, 



namely .R-parity violation ]fL p [16|, [ill, ll8| . Allowing for all possible gauge- and supersymmet- 



ric invariant terms within the MSSM particle content, -R-parity can be broken by bilinear 



or trilinear terms [17|. Trilinear models provide a huge number of free parameters, where 



different attempts tried to reduce them [19[. Spontaneous -R-parity violation can be under- 
stood as a violation of lepton number by the vacuum expectation value of some singlet field 
16l |20L |21| . Bilinear ]fi p can be interpreted as the low-energy limit of some spontaneous 



Jji p model, where the new singlet fields are all decoupled. The bilinear model has only six 
new 1/ip parameters and is thus more predictive in particular in view of collider physics once 
the parameters have been chosen to explain correctly neutrino data. For an early work on 



bilinear .R-parity violation in the context of a low-energy effective theory we refer to [22 



*A recent analysis combining the newest precision results from various experiments like KamLAND 
Super-K 0, SAGE 0, SNO @] or MINOS [g can be found in Ref. 0. 



The phenomenology of Jfi v SUSY has been studied extensively in the past, for reviews 



see [23|, |24j, |25J . Neutrino masses have been calculated with trilinear couplings [17J, |26[] and 
for pure bilinear models [27|, |28|, |29|, |30|, |31[ or both [32[ . Neutrino angles are not predicted 
in either schemes, but can be easily fitted to experimental data. In bilinear schemes the 
requirement to correctly explain neutrino data fixes all rft p couplings in sufficiently small 
intervals, so that in some specific final states of the decays of the LSP correlations with 
neutrino angles appear. This has been shown for a (bino-dominated) neutralino LSP in 



33l |34|, |35|, for charged scalar LSPs in [36 



37|, |38J, and for chargino, 



for sneutrino LSPs in 

gluino and squark LSPs in [37j. Such a tight connection between neutrino physics and LSP 
decays is lost to some extent in t he g eneral trilinear-plus-bilinear case. 

The /i-problem of the MSSM 39[ is induced by the mass term for the Higgs superfields 
fiHdH u . For phenomenological reasons the parameter \x must be of the order of the elec- 
troweak scale. However, if there is a larger scale in the theory, the natural value of \i should 
be around this large scale. A solution to this problem can be provided by the Next-to- 



Minimal SSM (NMSSM) [40l. |41| . where an additional singlet superfield is introduced. The 
vacuum expectation value of its scalar component generates an effective /i term as soon as 
electroweak symmetry is broken. The phenomenology has been discussed for example in 



42, 43 . Electroweak mass corrections were discussed in 44 , an on-shell scheme for neu 



tralino and chargino masses including the discussion of the two-body decays x°i ~~ >" xfW T 



at one-loop level was worked out in [45 



The /ii/SSM, which was first proposed in [46|, uses the same singlet superfield(s) not 
only to generate the /i term but in addition Dirac mass terms for the left-handed neutrinos. 
The cubic self-couplings of the singlets in the superpotential explicitly break lepton number 
and moreover also R p . Thus, Majorana neutrino masses are generated once electroweak 
symmetry is broken. One-loop corrections to the neutrino mass matrix have been calculated 
in 



47, 48 



Compared to the NMSSM the /xz/SSM with one right-handed neutrino superfield 
also contains six new IJip parameters, the same number as going from the MSSM to the 
bilinear ljl v model as will be discussed in detail later on. This implies that one also obtains 



correlations between neutrino mixing angles and ratios of LSP decay branching ratios [471 149 



50 



Moreover, one has unusual decay modes and signatures in the Higgs sector [471 . 

Although one expects correlations between neutrino mixing angles and final states of the 
LSP decays at tree level, it is a priori unclear what happens at one-loop level. In particular 
using one-loop on-shell masses and the one-loop mixing matrix together with the tree level 
decay widths results in case of the /i^SSM with one right-handed neutrino superfield in an 
unexpected behavior: The expected correlation does not necessarily show up for a singlino- 
like LSP. Therefore, we present in this paper a full one-loop calculation for the decays 
Xi —> l^W 11 and show how the tree level expectations evolve at full one- loop level. For 
this aim we extend the on-shell scheme presented in [45[ for the decays x°i xfW T within 
the MSSM and the NMSSM to the models considered here. Moreover, we briefly discuss 



the implementation of these decays in CNNDecays [5l|. Note that in case of the /iz/SSM 
with more than one right-handed neutrino superfield similar problems do not show up. The 
reason is, that neutrino physics can be fully explained at tree level and thus at the same 



3 



level of perturbation theory as the correlations with the neutralino branching ratios. 

This paper is organized as follows: We first define our models in Section [21 In Section [3] 
we present our on-shell scheme for the neutralino and chargino mass matrices in the consid- 
ered models. Afterwards we present numerical results for the neutrino masses and mixings 
comparing the on-shell scheme and the DR scheme in Section |H Next we discuss the decays 
x5 — > / ± W /=F and the correlations to neutrino physics at tree level in Section In Section 
we first discuss the complete NLO corrections to the widths and discuss the correlations at 
one-loop level. Finally we conclude in Section [7J Appendix A gives the implementation of 
the considered models in the public program CNNDecays. 



2. Model definition 

In this section we introduce the models by presenting the superpotential and the soft- 
breaking terms. Moreover we will give the neutralino and chargino mass matrices at tree 
level and comment on the need of one-loop corrections for the neutrino masses. 

2.1. Superpotential, soft terms and scalar potential 

We will start with the superpotential of bilinear i?-parity violation (BRpV) and the 
/ii/SSM with one right-handed neutrino superfield v c . Both model have in common the 
Yukawa part of the matter superfields coupled to the SU (2)-doublet Higgs fields 



where e a p is the complete antisymmetric SU(2) tensor with = 1 and i,j denote the three 
families. In case of BRpV the explicit [x and e^-terms are added and the total superpotential 
reads as 

(2) 



whereas in the /i^SSM the Higgs doublets are coupled to the right-handed neutrino super- 
field v c and in addition a term containing the neutrino Yukawa coupling Y v is added 



H — ki/V^ 
3 



(3) 



The first term in Equation (J2J) and the last two terms in Equation ([3]) explicitly break lepton 
number if one assigns lepton number to v c . If not the first term in Equation ([3]) explicitly 
breaks lepton number and one immediately sees that in both models one has the same 
number of Jjt p violating parameters. Note that for the phenomenology it does not matter 
if v c carries lepton number as it is broken explicitly by a least one interaction of this field. 
Both models have in common the soft SUSY breaking terms 



r 



MSSM-B^ 



soft 



m 



ij 2 icjc* 
i 3 



m 



2 H HTH a d + m\H a :E a u - i [M^B + M 2 W^W> + M 3 ~g^ + h.c] 

4 



(4) 



where K^f M B ^ contains all the usual soft terms of the MSSM but the 5 u -term and a 

soft F 

summation over 7 = 1, . . . , 3 and 7' = 1, . . . , 8 has to be performed. In addition we add in 
case of BRpV 



V. 



b$ P 



soft 



yMSSM-Bfi 
soft 



e a/3 



B^Ei - BvHSHl + h.c. 



(5) 



and in case of the /xz/SSM: 



T/ /^SSM _ T/ MSSM- J B M 2~c~c* 
V soft ~ V soft + m ^ c V V 



+ /i.e.] 



(6) 



In these expressions the notation for the soft trilinear couplings introduced in [52| is used. 
We will not present the scalar potential and the corresponding tadpole equations in this 
paper, but refer to 301 fo r BRpV and to 47| for the /i^SSM. Note the different definition 
of k in comparison to 47J. The determination of the vacuum structure through the scalar 
potential induces the following vacuum expectation values (VEVs): 



V2 



Vd, 



'X) 



7f 



7f 



7f 



(7) 



Again v c is only present in the //i/SSM, where its vacuum expectation value generates effec- 
tive /1- and ej-terms 



—Y l v 
y/2 ve 



(8) 



In this way v c plays the role of the singlet superfield S in the NMSSM, implying that /i, but 
also 6j is naturally of the order of the electroweak scale 43[. Moreover there are no additional 
terms with dimension of mass present in the superpotential. Beside the parameter 6j we 
define the so-called alignment parameters 



Aj = fAVi + V d 6i 



(9) 



which determine the neutrino mass matrix at tree level |47| and which give the relevant 
information on the tree level mass matrix in a basis invariant way 53] . 6j and Aj can then 
be used to fit the neutrino data. Concerning the scalar potential we add, that in case 
of BRpV we determine the soft-breaking parameters m 2 H , rn 2 H and Bi and in case of the 
/i^SSM m 2 Hd ,m 2 Hu ,mlc and T l v from the tadpole equations. 



2.2. Masses of neutralinos and charginos 
Using Weyl spinors defined by 



(ry 



W~,H. 



T 



w + , h: 



e c ,/i c ,r c 



(10) 



the Lagrangian density containing the chargino masses can be written in the form 

1 



C D 



((r) T M c ^ + + ^ + ) T (M C ) T ^-) + h.c. 



where the mass matrix of the charged fermions is given by 



The mass eigenstates are obtained by two unitary matrices U and V 

F+ = V it il>t and Fr=U lt ^ , 



( M 2 


T2 gVu 















Lyii,,. 


Ly» 2 ,]. 


V2*i v 


T2 gVl 


-ei 




>/5 r « Vd 




7f^ 2 






-±-Y 22 iu 

^2^1 Vd 


75^ 


\T2 9V * 


-C3 




V2*l Vd 


7^ 



i < j 



so that the mass eigenvalues can be calculated via 

M C 4ia. = Diag . . . 

m-+ < m -+ for 

and the Dirac spinor xt can be deduced from 

xt = ( i 

For the neutralinos we obtain in the basis H 
the mass matrices of the neutral fermions, which have the generic structure 

M n ~- 



M H m 
T 



m 



'ir. 



(12) 



(13) 
(14) 



(15) 



(16) 



(17) 



and enter the Lagrangian density C D — | ((ip° ) T Ai n ip°) +h.c. The submatrix Mh contains 
the 4(5) heavy states, which are the four usual neutralinos of the MSSM and in case of the 
/ii/SSM in addition the right-handed neutrino. The submatrix rh mixes the heavy states 
with the left-handed neutrinos and contains the i?-parity breaking parameters. In detail the 
elements are given by 



H 



I 


Mi 





-\g'v d 


\g'v u 


\ 







M 2 


\gvd 


-\gv u 







-\g'v d 


\gv d 







~T2 XVu 




\g'v u 


-\gv u 


-n 





~T2 Xvd 


V 








-^ Xl ' u 


~T2 Xvd 


y/2nv c J 



l v c is only present in the /Lti/SSM. 



6 



(-h'vi \gv x e l 



m 



2< 



g'v 2 \gv 2 e 2 e 2 ^ 
\-\g'^ \gvz e 3 e 3 fj 



(19) 



where in the /ii^SSM /1 and e, can be taken from Equation (JSJ) . In case of BRpV one has 
to omit the last column and the last row in the submatrix Mh and the last column in the 
matrix in. The mass eigenstates are obtained via 



F° = ATM 



(20) 



from the gauge eigenstates i/j®, where the unitary matrix M diagonalizes the mass matrix 
M n as 



M^dia. = Diag fmjo, 



N*M n M ] 



(21) 



where j = 7(8) in case of BRpV (/xi/SSM). The masses are ordered as m^o < m^o for % < j 
and we have chosen M such that all mass eigenvalues are positive. This implies that in 
general the matrix M is complex. Similar to 45[ we can do useful checks using negative 
mass eigenvalues in combination with a real mixing matrix N '. The 4-component spinor can 
be obtained via 

n 

F? 



"X,i 



(22) 



All the other masses of the scalar sectors can be taken from |30| for BRpV and from |47j for 
the /xz/SSM. The various couplings, which are needed for the one-loop calculations, can be 



found in the folder couplings within CNNDecays [51 



2.3. Neutrino masses 

Due to the smallness of the i?-parity breaking parameters an effective neutrino mass 
matrix in a seesaw approximation can be found 



m 



eff 



—m T M H 1 rh 



(23) 



where the matrix £ contains the small expansion parameters and can be taken from |30|, |4 
The basic calculation yields 



cff\ 
vvJij 



aAjAj 



4Det 



with a being 

^ssm _ m 7 (X 2 v d v u + (j,m R ) 
C 



4/iDetr SSM 



where we have used the abbreviations in Equations (jSJ) and (Q and 

g 2 M 1 + g' 2 M 2 , m R = s/2kv c , 



(24) 
(25) 

(26) 



where the singlet mass parameter rriR is only valid in the yU^SSM. The determinants of Mh 
are given by 



Det 



bijtp 



1 

-m^u,v d v u 



M x M 2 [x 2 
7 



(27) 



DeC SSM 



m 7 (A v + Am R /j,v d v u ) - M 1 M 2 ^(v d v u X + m R ji) 



f28) 



The projective form of this mass matrix implies that only one neutrino acquires a tree 
level mass, while the other two remain massless. We denote the matrix, which diagonalizes 
in the following M v . The correct magnitude of A can be easily estimated from the 



m 



off 



approximated formula for the neutrino mass m V3 at tree level 

a|A| 



m,. 



(29) 



Since one neutrino mass is not sufficient to explain the oscillation data, one-loop corrections 



in BRpV |30|, l3]J, [54( and in the /xz/SSM |47|, |48|, |55| have to be added to match the current 



bounds. We will therefore consider the complete mass matrices at the one-loop level using 
an on-shell scheme for calculating the masses. In addition we will show, that this also 
requires the complete NLO corrections to the decays T(x5 - ► I W^) to obtain the correct 
correlations between neutrino mixing angles and ratios of branching ratios. 



3. One- loop on-shell masses of neutralinos and charginos 



As we have shown pure tree level neutralino masses including the neutrinos cannot ex- 
plain the full neutrino bounds in BRpV and the /x^SSM with one right-handed neutrino v c 
only. We will therefore go to the next order in perturbation theory for the masses, which was 



illustrated using a DR scheme in BRpV |30|, |3l|, |54J and in the yUi^SSM |47|, |48_|, |55| m previous 
articles. However, we will later focus on the decay width T(xi ~> lfW T ) at one-loop level to 
investigate the relations between branching ratios and neutrino mixing angles on higher loop 
order and choose a more convenient on-shell scheme for our calculation. Before discussing 
this topic, we will first discuss the masses for neutralinos and charginos in an on-shell scheme 



as m 



45 



which allows to have finite corrections to the neutrino masses m}, L but also to 
the mixing matrix M 1L at one-loop level, which is crucial to explain the neutrino oscillation 
data. The approach is similar to 56|, |57|, whereas for the MSSM another procedure has 
been worked out in 58, 15911 . 



3.1. Common basic principles 

A special feature of the chargino/neutralino sector is, that the number of parameters is 
lower than the number of imposed on-shell conditions. This implies finite corrections to the 
tree level masses of neutralinos and charginos. We closely follow 45|, |56[ and start with the 
one- loop contributions to neutralino masses 5A4® and chargino masses 5M.f 



(6M$y = 6(Ar r M* dla .Ar) id 



Nni (SMl dia ) nl Mj+Ki (Ml dia ) nl 5Mj 



{5M% j = 8(U T Ml dia V) i 



SUr, 



III 



IJ 

V i:i + U n 



(30) 



(31) 



(tMl dia ) nl V lj + U n 
8 



(M 



c,dia) n [ SVlj 



with the diagonalized mass matrices M dia and their counterterms 8Ai dia '■ 

(Ml dia ) nl = 5 nl m 2 o, (8Ml dia ) nl = 8 nl 8m n (32) 
(Ml dm ) nl = S nl m^, (SM® dm ) nl = 5 nl 6m^ (33) 

The bare mass matrices of the neutralinos and charginos can now be expressed as the 
corrected on-shell mass matrix with the corrections in Equations (130]) and (}3~Tj) or via the 
tree level mass matrix expressed in physical parameters together with the renormalization 
constants of those: 

M° niC = M® c + 5Mt, c = M n , c + 8M n , c (34) 

Therefore the relations between tree level and one-loop mass matrices take the form: 

M% c = M n , c + 8M n , c - 8Ml c =: M n , c + AM n , c (35) 

In the following we will define the model-dependent physical parameters, write down the 
renormalization of the mass matrices 8Ai n ,c and identify the renormalization constants of 
the physical parameters, which should be fixed in the neutralino or chargino sector. Some 
physical parameters, namely Smw, 8m z and thus 8 cos 9w are fixed in the gauge boson sector 



and 8 tan f3 in the Higgs sector. In particular for 8 tan f3 we take the DR renormalization (60 
so that UV divergences in the masses and the considered process cancel 

= ^A (3Tr(y,y d t) - 3Tr(Y u Yj) + Tr(Y^t) - Tt(Y u YJ)) , (36) 

where A = - — 7 + ln(47r) parametrizes the UV divergence as in [45[ and Y v is only present 
in the /iz/SSM. Note that this choice maintains also the gauge invariance of masses and the 
considered decay widths. We will also need the renormalization constants for the sine and 
cosine of different angles, which are given by: 

5cos/3 . 2 p8 tan (3 8 sin (3 2 ^8 ^tan (3 
cos (3 tan/3 ' sin/3 tan/3 

3.2. Bilinear R-parity breaking 

The on-shell renormalization of neutralino and chargino masses in the MSSM and NMSSM 



was extensively discussed in (45| . In case of BRpV the renormalization of the physical pa- 
rameters is more challenging. The masses of the Z and W bosons are defined in the following 
form 

m 2 z = \{g 2 + g 2 ) [vl + vl + » m w = -/ Wd + v " + ^2 > ( 38 ) 

where not only the vacuum expectation values of the Higgs enter, but also the ones of the 
left-handed sneutrinos. As in the (N)MSSM the cos of the given Weinberg angle cos 9\y can 
be derived from mw = cos 9wmz- Beside the already known tan/3, we define 

v ■ 

tan f3i = — and 6j (39) 

v d 

9 



as additional parameters, which are present at tree level, such that the tree level mass matrix 
of the neutralinos is given by 



M H m 
T 



m 
/Mi 



M, 



H 





M 2 
sym. 



— rriz sin 9w cos (30 
mz cos Qw cos (30 




mz sin 6*^ sin (30 \ 
-mz cos #w sin (30 

/ 



(m T )i = (— mz sin^vK cos/3 tan $0 cos 6^ cos /3 tan/3j0 
where is a function of tan (3 and tan namely: 



e = e(/3,A) 



i 



1 + cos 2 /3 ^ tan 2 $ 



(40) 

(41) 
(42) 

(43) 



In this way we maintain the possibility to fix the renormalization constants of mz and cos 0\y 
in the gauge boson sector, whereas the corrections from i?-parity breaking are parameterized 
by tan/3j, and the parameter O, the latter depends on (3 and (3i. Defining in addition the 
lepton masses m^ = -^Y^Vd the chargino mass matrix is given by: 



( M 2 

\/2m w cos (30 
\/2mw cos (3 tan (3iO 
\/2mw cos (3 tan (3 2 

\\/2m w cos (3 tan /3 3 



\/2mw sin (30 
-ei 

-62 
-63 





tan Pima 
mn 

m 2 i 
m 3 i 



\ 

tan (3imi2 — tan (3im^ 

m 12 m Vi 

m 22 m 2 3 

m 32 m 33 J 



(44) 



jj, since the counterterms <5mij 



We discuss the general case of non-diagonal lepton masses m 
are needed to cancel non-diagonal contributions in the lepton mass matrix at one-loop level. 
However, the lepton Yukawa couplings Y\ and thus m^ can be chosen diagonal at tree level. 
Using now the relations cos (3° = cos (3 + 8 cos (3 and tan (3® = tan (3i + 5 tan (3i one can do a 
simple expansion in first order of 5 cos (3 and <5tan/3;: 



cos 2 /3°£.tan 2 /3° 
1 




l + cos 2 /3£.tan 2 A 
^^(cos 2 (3 tan (3iS tan (3i + cos /35 cos (3 tan 2 



f cos 2 (3 £\ tan 2 $ 

= + 50 (45) 
Therefore the counterterm of can be expressed as a function of 5 tan/3j and 5 cos /3: 

(56 = - cos (35 cos /30 3 tan 2 # - cos 2 /30 3 tan #5 tan A (46) 
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The variation of all non-zero entries of the tree level mass matrix yields: 

SM» = SM, = ^A-f" (47) 

MS = -6(m z ^B w ™m = f *2£ + ^ + ^ + ^ AC (48) 

\ m z singly cosp B y 

5^n 4 = 5 (m z sm6 w sm/39 = + , w + — -f + — KM* 4 49 

\ sin tV sin p By 

^M+i = _ 5 ( mz sin ^ cos tan ^.q) ( 50 ) 

<5m7 (Jsin^w <5cos/3 5 tan A 50 \ , , ■ 

+ • 7 + 7T + 7T + TV ^ 

sm Uw cos p tan pi By 

iAC = *M 2 = (51) 

M<? = ^cos^cos^e) = f *5» + ^ + *£5^ + ») A4^ 3 (52) 

\ costV cosp B y 

2 4 . . ( 5mz c^cos^iy 5 sin /3 50 \ 24 



ciMf = -5(m z cos^sin/30) = — + ^ + — -f + 77 M 4 (53) 

1 mz cos #14/ smp B 



5M 2 / +l = 5(m z cos 6 W cos (3 tan ft©) (54) 

Siriz t^cos^w 5 cos 8 5 tan A 5B\ . /9/U _,- 

+ + T + ~i IT + "7T" -W' 

cos^ cosp tanpj B J 

5M™ = -6n = 5 -^Mf (55) 
H 

SM*/ +t = 5e t = 5 -^Mf +l (56) 

whereas all the other variations 5.M* 2 = 5Ai^ = 5Ai^ 4+t = 5.M 44 = necessarily vanish. 
The corrections in the chargino mass matrix read: 

= 5M 2 = 6 -^Mf (57) 

6M? = V25(m w sin 80) = ( ^ + ^ + *%) (58) 

\ mw smp By 

5Mf = V25(m w cos 00) = ( ^ + ^ + ^ (59) 

\ m w cosp B J 

8Mf = 5(i= 6 -^Mf (60) 

6M? +i = 5{-Un0 k m ki ) = V - ^ (61) 

^ V tan/3 fc m ki J 
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SM 2 c +iA = V25(m w cos tan #6) (62) 
5mw S cos 8 S tan A SQ \ , 



mw cos/3 taxi Pi 6 / 

= = ^^,2 (63) 

= £ = ^^2+^+, (64) 

3 rriij 

and 5A^c' 2+l = vanishes. Similar to the (N)MSSM we fix S M± in the neutralino and SM 2 , S/j, 
in the chargino sector. However, we still have to find appropriate renormalization conditions 
for S tan Pi, Sei and Srriij. Whereas Srriij should be fixed in the lepton sector, so that one-loop 
contributions to leptonic two-point functions vanish, we have several possibilities for S tan Pi 
and 5ei. Whereas S tan /3, could for example be fixed in the Higgs sector together with tan/3, 
one could fix Se t with respect to an i?-parity violating decay. We will focus on stable lepton 
masses. Thus, we calculate Stan Pi and 8ci from the one-loop contributions in the chargino 
sector, which mix the well-known MSSM charginos with the leptons. In total we therefore 
impose the following conditions for the MSSM parameters 

AM 11 = AM 11 = AM 22 = , (65) 

which results in 

5M x = 5Mf\ 5M 2 = 8Mf\ 8fi = -8Mf 2 . (66) 

In a second step the conditions AAi 2+l ' 2+ ^ = for the renormalization constants of the 
lepton masses Srriij are imposed, resulting in 

5mii = 8Mf +i ' 2+j . (67) 

In a last step the i?-parity violating sector with 5 tan /3j and Sei is considered, which can be 
fixed by imposing AMl +h2 = AM. 2,2 ^ = 

1 1 X 

Stan Pi = - Y] e ijk ei rs Tim jr m ks (68) 

j,k,l,r,s 

with Ti = - tan 0kSm ki - 5Mf' 2+i (69) 

k 

5e l = 5Mt w , (70) 

where €ij k is the Levi-Civita symbol and det(m) is the determinant of the lepton masses. In 
case of vanishing non-diagonal lepton masses at tree level m^- = for i ^ j this simplifies to 
Stan Pi = — Tj. Another possibility is to fix 5 taxi Pi from AAi 2 ^' 1 . However, this induces 
a dependence on the renormalization constants of the gauge sector, whereas in the described 
case the neutrino and lepton sector are "decoupled" from those. 
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For all the remaining entries of the neutralino and chargino mass matrices shifts AA/f„ iC 
have to be taken into account. Due to the non-vanishing entries AAil +l ' 1 also the lepton 
masses differ between tree and one-loop level. However, we will show that this difference is 
tiny, for reasonable neutrino masses far below the experimental uncertainties. Note that the 
Yukawa couplings of the leptons at tree level of course have to be adopted, so that the tree 
level lepton masses fit the experimental known values. 



3.3. fiisSSM 

Combining the already known procedures from the NMSSM and BRpV we define 

Vi n v c 

— , tan p c = — 
Vd v u 



tan f3i 



1 



m R = ~^ KV a 



In this notation the tree level mass matrix of the neutralinos is given by 
M n -- 



M, 



H 



[m 



M H m 

* T 

m 

( Mi —mz sin 9w cos (3Q mz sin 9w sin f3Q 

M 2 mz cos^iy cos (3 6 — mz cos Ow sin (3Q 

-n 

sym. 

V 

— m z sin9 w cos/3tan/3j0 mz cos 9w cos /Stan faQ e. 








tan /3 C 
tan ft dj- fi 
tan p tan p c 



m R ) 



tan p c 



(71) 
(72) 

(73) 

(74) 
(75) 



where G is defined in Equation (143]) and the chargino mass matrix has the same form as 
in BRpV, see Equation (1441) . Beside the variations already present in the bilinear model as 
shown in Equations (f47|) - (l56|) . where the indices 4 + i have to be shifted to 5 + i, 5M. n has 
the following additional entries: 



5M 



85 



5 



15 



5M 



5Mf 



tan j3, 
( tan Pi€i - ii 
\tan (3 tan (3 C 

8m R 



5fi 



S tan f3 c 
tan j3 c 
5 tan P 



Mf 

5 tan P c 



(76) 



tan (3 tan (3 C 



6 tan j3i€i + tan ftiSei — 5fi 
tan f3i€i — ji 



M 



15 



5M 



5,5+i 



m R 
= 5 



tan /3 C 



5q 



S tan f3 c 
tan /3 C 



M 



55+i 



(77) 
(78) 

(79) 



whereas all the other variations 5M^ = 5M^ 5 = 5Mf = 5M%* = 5M^ 5+i = 5M^ = 
necessarily vanish. The corrections in the chargino mass matrix are the same as in the 
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bilinear case, presented in Equations ( |57|) - (|64|) . Similar to the MSSM, we fix 5 Mi in the 
neutralino and <5M2,5/i in the chargino sector. Similar to the NMSSM, we fix 5tan/3 c and 
5m a in the neutralino sector and similar to the BRpV case <5tan/3j and 5ei are determined 
in the chargino sector. However, we will summarize these results once again for the /wSSM. 
We start with the conditions for the non-i?-parity breaking variables 



AM 



11 



AM 



22 



AM 



ii 



AM 



AM 



55 







30) 



which results in: 



5M 1 

S tan f3 c 



5Mf\ 
tan 2 /3 C 



5M 2 
SMf 5 - 



5Mf\ 
1 



8fi 



-5M 



©34 



tan f3 c 
-- AM 2 : 2+l 



5Mf 2 
Sm R ■■ 



5M 



©55 



(82) 



In a second step the conditions AM 2+4,2 = AMl' 2+t = AM 2+t ' 2+j = for the renormaliza- 
tion constants of the lepton masses dm^, 5tan/3j and Se^ are imposed, resulting in: 

be t = 5Mf +l ' 2 (83) 

(84) 



Srriij = 5Mf +i ' 2+j , 



5 tan 0i 



det(m) 2 
with Tj = 



k^-lrs^- l^^jr^'ks 



j,k,l,r,s 

tan [3kSm 



hi 



5M 



©2,2+i 



55) 



3.4- Definition of one-loop on-shell masses 

With the procedure introduced in the last sections we can calculate one-loop on-shell 
neutralino and chargino masses for BRpV and the /iz/SSM, namely by combining the full one- 
loop corrections 5M® C with the counterterms 5M n>c obtained according to Equation (155]) . 
This results in the one-loop mass matrix M® c , whose diagonalizations lead to one-loop neu- 
tralino mlo and chargino masses ml± and mixing matrices at the one- loop level Af lL , U 1L , V 1L 
Note that these masses are UV and IR finite as well as gauge independent, if one takes into 



account the gauge independent renormalization of the mixing matrices as presented in |45 
for the Equations (130]) and ( 131]) . Similar to the case of DR masses an effective neutrino mass 
matrix of the form 



[m 



eflf 



)ij = aAiAj + b(eiAj + ejAi) 



cei€j 



(86) 



is obtained at one-loop level |30|, [31 



3.5. Calculation of the neutrino parameters 

With the one-loop on-shell masses m^ L of the neutrinos and the mixing matrix M 1L one 
can calculate the relevant parameters to be compared with experimental neutrino data: 



Am 



atm 



tan 2 9 atm 



\f iL 

- /v 3,6 



Kf) 2 > 



\f 1L 
- /v 3,7 



tan 2 # 



57) 



sol 



Kf 1L 

■ /v l,5 



U, 



c3 



lip 

3,5 



The formulas are valid for BRpV, in case of the /xz/SSM one has to replace Aflj — > A/"^ 7 L „ 



'i,j+V 
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Parameter 


Best fit 


2(7 


Am^lO" 5 eV 2 ] 
lAmlilflO- 3 eV 2 ] 
tan 2 6\2 = tan 2 9 so i 
tan 2 6*23 = tan 2 atm 
tan 2 #13 = tan 2 6 R 


4^+0.09 

U.^±dO_ 03i 

1 04 +0 - 28 
n nin+ 009 


7.24 - 7.99 
2.28 - 2.64 
0.39 - 0.54 
0.69 - 1.56 
< 0.028 



Table 1: Current bounds on neutrino data taken from [1 



4. One-loop on-shell neutrino masses 

In this section we show the results for the one-loop on-shell neutrino masses in BRpV 
and the /xi^SSM, which we obtain using the definitions of the last sections. The two models 
and the procedure explained before are also added to CNNDecays [1H, so that the reader is 
able to reproduce all the presented results. 

4-1- Parameter choice 

In the choice of our parameters we closely follow the strategy presented in 47]. In Tabled] 



we show the current bounds on neutrino data taken from 10|, which will be used for our 
later analysis. Moreover we note for the specific models: 

• BRpV: We use the low-energy parameter sets of the well-known "Snowmass Points 



and Slopes" (SPS) |6l|, |62| as benchmark scenarios. In the discussion of the decay 
widths we refer to two particular scenarios, which are defined in Table [2j SPS 2 
was changed such, that a higgsino-like lightest neutralino is present by setting Mi = 
M 2 = 600 GeV. The additional parameters e^Vi or respectively e i5 Aj are fixed by the 
neutrino constraints, which can be found in Table [TJ if not mentioned to be chosen 
otherwise. The corresponding Bi are determined from the tadpole equations. 

/zz/SSM: We refer to the low-energy parameter sets of the NMSSM benchmark scenarios 



c 



[42l . |63| named mSUGRA i or GMSB j based on the soft SUSY breaking mechanism. 
Details can also be found within jiij. We add for completeness, that the VEV v, 
is fixed from /i and A using Equation ([8]). In the discussion of the decay widths we 
use four scenarios, which are defined in Table |2] and which are chosen such, that the 
lightest neutralino is either bino, wino, higgsino or singling. Concerning perturbativity 



and charge or color breaking minima we refer to |47l . |64| and note that our choices 
of A, k and T\,T K respect the given constraints. The additional parameters Vi,Y* or 
respectively Aj are fixed by the neutrino constraints, which can be found in Table (TJ 
if not mentioned to be chosen otherwise. The corresponding T* are determined from 
the tadpole equations. 
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SPS 2' 


SPS 3 


mSUGRA 1 


mSUGRA 3' 


mSUGRA 4' 


GMSB 5 


Changes 


Mi = M 2 
= 600 GeV 


None 


None 


Mi o M 2 


/c ~ 0.1 
-> 0.4 


None 


Xi ■ m 

1 L 

m 
C 


371.05 
370.68 
H 


165.03 
164.76 
B 


210.79 
210.61 
B 


208.78 
208.55 
W 


196.82 
199.51 


203.30 
203.30 

S 


xl ■ m 

1 T, 

m 
C 


395.76 
395.67 
H 


295.36 
295.73 
W 


387.18 
387.10 
W 


391.11 
390.82 
B 


205.60 
205.24 
H 


496.87 
496.81 
B 


xl ■ m 

1 L 

m 
C 


600.00 
599.87 
B 


504.42 
504.76 
H 


971.11 
971.75 
H 


941.92 

r\ A r\ r\r\ 

940.00 
H 


327.26 
326.83 
5 


899.60 
899.98 
W 


xl ■ m 

1L 

m 
C 


624.71 
624.42 

w 


521.86 
520.49 
H 


976.52 
975.14 
H 


942.49 
943.05 
H 


330.44 
330.95 
B 


1377.67 

1 o 'i'i a r 

1377.45 
H 


xl ■ rn 
C 






2101.57 
2101.57 

5 


1421.70 
1421.66 

S 


608.43 
607.65 
W 


1383.97 
1383.04 
H 


Xi - m 
C 


377.77 
378.22 
H± 


295.07 
295.58 

w± 


387.16 
387.23 


208.83 
208.77 


201.36 
201.75 
H± 


899.59 
900.14 


xt - m 

C 


619.84 
619.13 
W ± 


521.90 
521.03 


977.07 
976.69 
H± 


946.04 
945.71 
H± 


608.40 
607.81 
W ± 


1384.24 
1383.51 



Table 2: 2 scenarios for BRpV and 4 scenarios for the /i^SSM to allow for lightest neutralinos with different 
particle characters C where m and m 1L represent the tree level and one-loop masses; SPS 2, mSUGRA 3 
and 4 were adopted, so that also a wino- and higgsino-like lightest neutralino is accessible. 



16 



4-2. Absolute neutrino masses and mass differences 

First we illustrate the behavior of the absolute neutrino masses and point out, that the 
on-shell renormalization allows a similar parameter dependence as the DR masses, which 
were defined in 3^ • 



> io 1 
■B io° 

io- 1 

IO" 3 
1(T 4 

io- 5 
io- 6 

IO -7 
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a) 



\ 
I 
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k17|A| 



Lf 10 1 r 



IO' 2 
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Figure 1: Three on-shell neutrino masses m v in eV as a function of |e] /|A| for a scenario in a) the /i^SSM 

> io 2 1 



based on mSUGRA 1; b) BRpV based on SPS 3 
> io 




Its 10 1 = 



io" = 



IO" 1 r 
IO" 2 r 
IO" 3 r 

10 -4L_L 



b) 



kf/|A| 



10" 



Figure 2: Three on-shell neutrino masses m\, L in eV as a function of |e| 2 /|A| using the sign-condition 
defined in Equation ([55)1 for a scenario in a) ^SSM based on mSUGRA 1; b) BRpV based on SPS 3. 



Figure [T] shows the neutrino masses of the three left-handed neutrinos as a function 
of |e| 2 /|A|, where we have chosen e\ = e 2 = e-s and Ai = A 2 = A 3 and a fixed value of 
|A| = 0.235 GeV 2 . In case of the /iz/SSM the scenario is based on mSUGRA 1 and in case 
of BRpV on SPS 3. In Figure [2] we set Ai = — A 2 = A 3 , so that the sign-condition 

— -— — < 89 

e 3 As 

is fulfilled. This allows for a simpler fit to the solar angle, since it helps to decouple the 
atmospheric and the solar problem by reducing the contributions from the b-term in the 
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effective neutrino mass matrix at one-loop level in Equation (186|) 30|, |3l|. We will therefore 
make use of the sign-condition in the following. Both models show a similar behavior. 

The absolute value of |e| determines the neutrino masses m]^ and m]^, which are gen- 
erated at one-loop level, whereas the tree level neutrino mass m]^ is set constant by |A| 
for small |e| and is only affected by the one- loop corrections for large values of |e|. Note 
that for the explanation of the neutrino data the individual q and Aj have to be chosen 
differently. In case of the scenarios without sign-condition a level-crossing as in Figure Q] a) 
can take place, which existence is dependent on the parameter point in the /iz/SSM as well 
as in BRpV. It corresponds to a sign-flip between m]^ and m]^. 

The Yukawa couplings Yi of the leptons at tree level have to be adopted such, that the 
tree level lepton masses coincide with the experimental values. Suppose we have fitted the 
tree level lepton masses mi to the experimental values, we can define the relative one-loop 
corrections by 



Si 



m 



IL 



mi 



rni 



(90) 



For the /x^SSM scenarios already presented in Figures [T] a) and [2] a) with respect to the 
neutrino masses, Figure [3] shows the relative correction 5i for the lepton masses at one- 
loop level. For reasonable neutrino masses the corrections to the lepton masses are small 
5i < 10~ 10 , even for the electron below the experimental uncertainties. Note that the shown 
dips correspond to a sign change, since we show the absolute value of the correction. 



t- m-ii 
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Figure 3: Corrections 5i denned in Equation (|M| for the three lepton masses (r (black, solid), /i (red, 
dashed), e (blue, dot-dashed)) as a function of |e| 2 /|A| for the /ii^SSM as in Figure [1] a), in detail a) without 
sign-condition ([55)1 ; b) with sign-condition 



In the following we want to point out the difference between the DR masses in [30J and 
the on-shell masses as defined in this article. For BRpV based on SPS 3 we present both 
schemes in Figure HI 

Although there exists a difference in the lightest neutrino mass m^, the mass differences 
lS.m\ tm and Am^ defined in Equation (1H7|) are very similar, since those are given by the 
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Figure 4: Mass differences Am„ tm (black, solid) and Am^, (red, dashed) as a function of |e| 2 /|A|for BRpV 
based on SPS 3 using the a) on-shell masses defined in Section [3] b) DR masses as defined in 



absolute values of m]^ and m]^. Since the on-shell procedure tends to reduce the impact of 
the one-loop contributions, the largest mass m]^ is in some scenarios only for larger values of 
|e| affected by the one- loop contributions when comparing the on-shell with the DR masses. 
However, for the neutrinos the two renormalization prescriptions are very similar, whereas 
for the other neutralinos and charginos and in particular for the leptons the corrections are 
(much) smaller in case of on-shell masses compared to the DR masses. In the scenarios above 
5i for DR lepton masses is between 0.3 and 2%. The one-loop corrections to the neutralino 
and chargino masses are included in Table |2] and are typically in the per-mil range, whereas 
in case of DR masses the corrections are in the order of a few per-cent 44 . 



4-3. Relation between A, e and the neutrino mass differences/the mixing angles 

In this section we discuss the relations between the neutrino mass differences, the mixing 
angles and the alignment parameters e and A. We refer to 10] for the current neutrino data, 
which we presented in Table [TJ For all figures presented in this section we fit e and A such, 
that the neutrino data bounds are fulfilled except for the angle shown in the figure. 

Similar to 30j we use A at tree level to explain the atmospheric mass difference and the 
atmospheric mixing angle, whereas e respectively e at one-loop level is used for the solar 
mass difference and the solar mixing angle. The latter correlation to e is more distinct than 
the one to e, where a prerotation with the matrix M v diagonalizing the tree level neutrino 
mass matrix rrfjl given in Equation (|24p according to e = M^e was performed. The mixing 
matrix M v is in the considered scenarios exactly given by 



65 



AC 



|A| 



\ |A| 



AiA 2 



VAl+Apl 
Aj 

v/Af+Af 
A2 

|A| 



A1A3 



\ 



V^Af+Apl 
A 2 

y/Al+Al 

A3 

|A| 



(91) 
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We will first comment on the atmospheric and solar mass differences Am^ m and Am 2 sol . 
Whereas the atmospheric mass difference is determined by |A|, the solar mass difference 
correlates with |e|. Therein the absolute value of A can be estimated using Equation (12*?]) . 




Figure 5: Correlation between the alignment parameters and the neutrino mixing angles for the /iz/SSM based 
on mSUGRA 1 (black), 3' (red, dashed), 4' (blue, dot-dashed) and GMSB 5 (brown, dotted) presented in 
Table [5J in detail a) tan 2 d a tm as a function of A2/A3; b) tan 2 8 so i as a function of Ixjl^- 

Using the definitions of Equation (1551) for the neutrino mixing angles we will now focus 
on the correlation between the alignment parameters A, and e, and the neutrino mixing 
angles, which is shown in Figure |5j The ratio A2/A3 defines the atmospheric angle and 61/^2 



fits the solar angle. The reactor angle is determined by A 1 /y / A|~+~Aj 30] which however 
can receive sizable loop corrections. 



5. Tree level decay widths of x° — > Z ± W r=F 

After having discussed the finite mass corrections between tree level and one-loop level, 
which allow to explain the full neutrino spectrum, we will now turn to the decays Xi ~~ >* 
I in order to study the correlations between neutrino mixing angles and the different 
leptonic final states. The partial widths for the decays x°j ~~ xfW^ including the ^decays 
Xi —> lfW T are obtained from the following interaction Lagrangian 

£ = ^T7 M {PlOuj + P R Rlj ) x°jW- + h.c. . (92) 
The couplings are given by 



Om = -gM; 2 U n - -Lg Lf*p a + g ^, 2+ ^ 4{5)+fe j (93) 
Omj^-gV^+^gV^ , (94) 
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where in Equation (193|) the index is 4 (5) in case of BRpV (/xi/SSM). The widths have the 
form 

r° = Y^^irnlrnlrn^-^lMTl 2 , (95) 



pol 



where m 8 (m Q ) is the mass of the mother (daughter) particle and Mt is the tree level matrix 
element. Explicitly they are given by 



K[m 2 -, 771?, 771? 



r° {x°j -»• xt w ) - 167rm s v "^"'i' ' '"w 

X ((|0 Wj | 2 +|O raj | 2 )/(771?,77lf,771^) - 6Re(0 Z//j Ora j )771 J ?71 i 



k(t71 2 , m 2 ,m 2 w ) 



x ((lOiifc^+IOffifcl 2 )/^ 2 ,^,^) _ 6Re(0 Lifc 0^ fc )77ii?7i fc 



using the functions 



, 1. . (x-y) 2 
f\ x i V, z) =- (x + y) - z + — — — 

k(x, t/, z) =x 2 + y 2 + z 2 — 2xy — 2xz — 2yz 



(96) 
(97) 

(98) 
(99) 



Coupling Xi — If ~ W ± 

The tree level couplings Oui in Equation (19U1) and Orh in Equation f[9"4"|) for the special 
case of the lightest neutralino x\ coupling to a lepton l± can be approximated in terms of the 
alignment parameters Aj and 6j, which allows to understand the correlations to the neutrino 
mi xing angles. Although we have shown these approximated couplings in great detail already 



in 



471 ] for the /iz/SSM, their knowledge is of great importance. The right-handed couplings 
Ofjn can be approximated in BRpV and the /ii/SSM by: 



Orh 



g(Yi)iiV d 
2Det+ 



gv u N 12 - M 2 N M g(2fi 2 + g 2 v 2 u )N 12 - g 2 {v d ^ + M 2 v u )N u ' 



fi 2/iDet + 
Therein the determinant of the heavy states in the chargino mass matrix appears 

Det+ = M 2 /i - ^g 2 v d v u . 



(100) 



(101) 



In contrast the left-handed couplings differ in both models and are given by the following 
expressions: 



n bF/L p ^9^i 
U1 ~V2 



g'M 2 fi 



2Det 



bljtp 
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fiDet + 
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N 12 

Vdm^ 



4Det 



b$ p 
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(102) 
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g'M 2 

2Det ^SSM 



(v d v u X 2 + m R jj)N_ 
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g gM 



Det+ 2Det^ SSM 



'12 
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8/xDet 



m. 



7 



{X 2 v d (v 2 d + <) + 2m R ^v u ) N 13 



V2/iDet + 8/xDet^ SSM 

777/ 

)et Q 
Am 



(103) 



4V2D<" SSM 

The couplings include the abbreviations of Equations (JHJ) , (JHJ) and f l26|) and the determinants 
shown in Equations (I27|) and f[28l) . Neglecting the right-handed couplings, which are gener- 
ally smaller than the left-handed ones, it is obvious that the lightest neutralino Xi couples 
to lfW T proportional to Aj without any dependence on the parameters. This fact is 
in case of the considered BRpV and the /iz/SSM independent of the nature of the lightest 
neutralino Xi- However, it is important to emphasize that all previous formulas are tree 
level results. A priori it is not clear, what happens if one uses the mixing matrices J\f, U and 
V for masses at one-loop level or performs a complete NLO calculation of the decay width. 

6. NLO decay widths of x° — > Z ± VF =F and correlations to neutrino mixing angles 

In this section we focus on the NLO decay widths T 1 (xi l + W~). The calculation 



procedure was introduced in |45[ and can be adopted for BRpV and the /ii/SSM without 



any changes. The corresponding formulas have been implemented in the public program 



CNNDecays [51| . All technical questions concerning the used renormalization scheme like the 
cancellation of UV and IR divergences or gauge invariance have been discussed in jiij for 
the (N)MSSM and preserve their validity for the i?-parity violating extensions. 

6.1. Absolute size of the corrections 

In Figure E] we present the NLO decay width for a mSUGRA 4 based scenario, where we 
varied k — 0.1 — 0.5 in order to have a singlino- and a higgsino-like lightest neutralino Xi 
in comparison. For k < 0.1 the decay is near the kinematical threshold and for k > 0.5 a 
Landau pole appears. We show the relative correction S(xi — > l + W~), which is defined by: 

«x? -» tw~) = r'(x;^+n--)-r°(x?^+^-) 
* ' r°oa-w+w-) y ' 

The size and the sign of the corrections is strongly dependent on the size and the sign 
of the i?-parity violating parameters and Y* or Aj and q. In Figure [6] we set A = 
(0.31,5.21,2.02) • 10~ 2 GeV 2 and e = (7.49,9.61,-6.57) ■ 10" 3 GeV. Since the solar mass 
difference and mixing angle are only well defined at one-loop level, the corrections to the 
decay Xi ~ > g + W~ are sizable and can be of the order of the tree level decay width or larger. 
The corrections to Xi ~^ H + W~ or t + W~ are generally smaller, but remain of order 10%. 
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Figure 6: Neutralino masses, particle content of the lightest neutralino and decay widths — > l + W ) as 
a function of « within the ^i/SSM based on the mSUGRA 4 scenario and A = (0.31, 5.21, 2.02) • 10~ 2 GeV 2 
and e = (7.49,9.61, —6.57) • 10 -3 GeV: a) Neutralino masses m^o (m~o (red, solid), mlg (black, dashed)); 

Xi Xl X-2,3,4 

b) Particle character \Afl^\ 2 {v c (red, solid), H u + Hd (blue, dashed), B (black, dot-dashed), W (orange, dot- 
dashed)); c) NLO decay width T 1 (e (blue, dot-dashed), /x (red, dashed), r (black)); d) Relative correction 
5 defined in Equation (|104[) (e (blue, dot-dashed), fi (red, dashed), r (black)). 

Correlation between decay widths and neutrino mixing angles 

Having discussed the absolute size of the NLO contributions to the decay widths of the 
two body decays Xx ~^ l + W~, we will now focus on the correlations between the branching 
ratios and the neutrino mixing angles. We have shown in Section that the tree level cou- 
plings can be approximated in terms of the alignment parameters A and e. Since the align- 
ment parameters are connected to the neutrino mixing angles as explained in Section I4.3[ 
one can expect the following tree level relation: 



r°ix^mv-)*{o^) - U 1 s>ton '-" (105) 



It is a priori not clear to which extent this relation also holds at one-loop level. We have 
found that using one-loop corrected DR masses together with the corresponding one-loop 
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Figure 7: Ratios T(x? -> /i + W _ )/T(x$ -> r+W - ) with r = T° in black, T = T 1 in red, dashed, T = T 
with 7V 1L in blue, dot-dashed, T = T° with M 1L ,U 1L and in brown, dotted as a function of A2 / A3 for 
the nvSSM scenarios of Tabled a) mSUGRA 1; b) mSUGRA 3'; c) mSUGRA 4'; d) GMSB 5. 

corrected mixing matrix M lL in the tree level decay width can spoil this relation. However, 
performing the full one-loop corrections including the real corrections from photon emission 
in the considered on-shell scheme retains the predictions at tree level, so that 



^ r+w _ } oc ^-j ~ tan 6 atm 



(106) 



In Figure [7] we highlight this feature using the scenarios of Table [2] for the /ii/SSM. We 
do not only show the tree level and one-loop relations, but also the relations obtained by 
using the one-loop mixing matrices Af 1L , U 1L and V lL within the tree level decay width 
r°(Xi — y l + W~), see Equation (196]) . The ratios are shown as a function of A2/A3, which is 
equivalent to varying tan 2 8 atm as can be seen in Figure |5j We vary the other parameters 
such, that the neutrino data are within their 2a bounds as given in Table [TJ This is the 
cause of the observed widths in case of the mSUGRA 4' scenario. 

As already known from BRpV the approximated one-loop contributions using T together 
with the one-loop corrected mixing matrices show the expected behavior if the lightest 
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neutralino is either mainly a bino, wino or a pure higgsino. However, in case that the singlino- 
component gets sizable or dominant, these approximations fail and the complete one-loop 
decay width T 1 has to be calculated to obtain a reliable result. In case of mSUGRA 4', 
which is shown in Figure [7J the higgsino has a purity of 91.7%, but a singlino component of 
6.9% is left, so that also in this case the full one-loop decay width is advisable. 

The reason for the misleading ratios in case of a singlino-like neutralino is a cancellation 
between the second and the last terms of the left-handed coupling Olii in Equation f )93|) 
at tree level. Replacing the tree level mixing matrix M by the one-loop mixing matrix 
M 1L spoils this cancellation, so that unreasonable results for the decay widths and also the 
ratios are obtained. However, the very same effect is restored once the complete one-loop 
corrections are included. 

For completeness let us briefly discuss the situation within BRpV. As examples we take 
scenarios where the lightest neutralino is either mainly a bino, scenario SPS 3, or mainly 
a higgsino, scenario SPS 2'. The ratios of decay widths using the various approximations 
are shown in Figure [SJ We observe that NLO corrections are somewhat more important for 
higgsino-like neutralinos: Here the combination of tree level decay widths with NLO masses 
and mixing matrices can lead to correlations which are off by a factor 2 compared the case 
of the complete NLO calculation. Moreover, the complete NLO corrections can shift the 
correlations up to 30 per-cent compared to the tree level result. 




a 2 /a 3 



10"' 10° 10' 10"' 10° 10' 

Figure 8: Ratios T(x? -> M + W~)/r(x? -> t+W-) with V = T° in black, T = T 1 in red, dashed, T = T° 
with Af 1L in blue, dot-dashed, r = T° with N 1L ', U 1L and V 1L in brown, dotted as a function of A 2 /A 3 for 
the BRpV scenarios of Tabled] a) SPS 2'; b) SPS 3. 
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7. Conclusion 

In ^models, where the sneutrinos obtain a VEV, a mixing between charged leptons and 
charginos as well as a mixing between neutrinos and neutralinos occur at tree level. The 
latter can be used to obtain at least one neutrino mass at tree level. The other neutrino 
masses can be generated at the one-loop level as it happens e.g. in bilinear ]fl p or in the 
/ii/SSM with one generation of right-handed neutrinos. In previous work a DR scheme had 
been used for this. In this paper we worked out an on-shell scheme and compared both 
schemes. The challenge here is, that there are less parameters than masses implying finite 
shifts at the one-loop level to the masses. With respect to the masses of neutrinos and 
charged leptons we have found only small differences well below the current and near future 
experimental accuracies. 

In both models one expects correlations between ratios of the decay widths or equivalently 
the branching ratios for the decays Xi ~~ >* l^W T and the atmospheric mixing angle which 
can be shown easily at tree level using approximate formulas for the mixing matrices. We 
have calculated these decays including the complete NLO order corrections in the on-shell 
scheme and compared the results to the pure tree level as well as to the case where one uses 
tree level formulas for the decay widths but one-loop corrected masses and mixing matrices. 

We have found that corrections to the channels containing either a fi or a r lepton are 
typically of order ten per-cent whereas the corrections for the e-channel can be of the order 
of the tree level decay width. The reason for the latter finding is that the solar mixing 
angle is not well defined at tree level but only at the one-loop level. For the correlations 
between the decay widths and the atmospheric mixing angle we find that in case of gaugino- 
like neutralino in all approximations nearly the same correlations. However, already for 
doublet-higgsino-like neutralino the correlations can be off by a factor two if one combines 
tree level widths with one-loop corrected masses and mixing angles. In case of a neutralino 
which has either a sizable singlino-component or which is dominantly singlino-like this is 
even more pronounced and the correlations get completely wrong in this approximation. 
Here one has to use the complete NLO corrections to obtain a reliable prediction. Last but 
not least we note that both models have been implemented in the publicly available program 
CNNDecays for the numerical evaluation. 
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Appendix A. New blocks in CNNDecays 



Both considered models, BRpV and the //i/SSM with one right-handed neutrino super- 
field, are now supported by CNNDecays, which was published in 45[ and can be found online 
0. We provide input files for the MSSM, the NMSSM, BRpV and the //i/SSM with one 
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right-handed neutrino superfield in the folder /examples, where in addition the correspond- 
ing output files can be found. In BRpV 6j and Vi can be chosen beside the input parameters 
of the MSSM: 



Block MODSEL 
1 
4 1 

BLOCK RVKAPPAIN 

1 1.45660382E-02 

2 9.01765562E-03 

3 -3.16131217E-03 
BLOCK RVSNVEVIN 

1 -8.68089903E-04 

2 -4.50162251E-04 

3 4.19592513E-04 



# Select model 

# bilinear model 

# RPviolation 



# kappa_l 

# kappa_2 

# kappa_3 



eps_l 
eps_2 
eps_3 



# 
# 
# 



_L_1 
.L_2 
.L_3 



The input file of the /iz/SSM follows the convention of the NMSSM, implying that the 
VEV of the right-handed sneutrino v c is calculated from the effective n = ^Af c . In addition 
k and the soft-breaking terms A\ and A K are input. The -R-parity violating parameters 
should be provided via the variables and vf 



Block MODSEL 

3 6 

4 1 

BLOCK EXTPAR 

61 1.00000000E-01 

62 1.08485437E-01 

63 -9.59966990E+02 

64 -1.58051889E+00 

65 9.68523016E+02 

73 5.71920838E-06 

74 6.20206893E-06 

75 -6.20206893E-06 
BLOCK RVSNVEVIN 

1 -1.35445088E-03 

2 -1.27271724E-03 

3 1.71364110E-03 



# Select model 

# munuSSM 

# RPviolation 



# lambda 

# kappa /2 

# T .lambda/ lambda 

# T_kappa/kappa = 

# mu_eff 

# Y.nu-1 

# Y.nu.2 

# Y.nu.3 

# v.L.l 

# v.L.2 

# v.L.3 



= AJambda 
A_kappa 



The corresponding soft-breaking terms of Y*, namely T^, are calculated from the tadpole 
equations in addition to m 2 H , rn 2 Hu and m^ 2 . The output does not only contain the im- 
parity violating decays of the lightest neutralino Xi ~^ , but in addition the relevant 
parameters for neutrino physics: 



DECAYTREE 
# 



1000022 

BR 

5.00974421E-01 
4.94172793E-01 
4.85278632E-03 

DECAY 1000022 

# BR 

4.99329442E-01 
4.95513455E-01 
5.15710289E-03 



2.03988581E-13 # chiOl 



NDA ID1 

2 15 

2 13 

2 11 
2.17956335E-13 

NDA ID1 

2 15 

2 13 

2 11 



ID2 

-24 # BR( chiOl 

-24 # BR( chiOl 

-24 # BR( chiOl 

# chiOl on NLO 
ID2 

-24 # BR( chiOl 

-24 # BR( chiOl 

-24 # BR( chiOl 



-> tau+ W-) 
-> mil + W-) 
-> e+ W-) 



-> tau+ W-) 
-> ran + W-) 
-> e+ W-) 



27 



Block SPhenoRP # additional RP parameters 



1 


5 


53918495E-02 


# 


eps_l 




2 


6 


00684651E-02 


# 


eps_2 




3 


-6 


.00684651E-02 


# 


eps_3 




4 


2 


10857140E-02 


# 


Lambda_l = 


v.d 


5 


2 


12780743E-01 


# 


Lambda_2 = 


v.d 


6 


2 


14264164E-01 


# 


Lambda_3 = 


v.d 


7 


2 


24691592E-03 


# 


m~2_atm [cV 


"2] 


8 


7 


74812242E-05 


# 


m~2_sol [eV 


-2] 


9 


9 


27311984E-01 


# 


tan "2 theta. 


.atm 


10 


■1 


08412270E-01 


# 


tan "2 theta 


.sol 


11 


3 


54098593E-03 


# 


U.e3"2 




15 


2 


40631533E+01 


# 


v.d 




16 


2 


40631533E+02 


# 


v_u 




17 


-1 


.35445088E-03 


# 


v_Ll 




18 


-1 


.27271724E-03 


# 


v_L2 




19 


1 


71364110E-03 


# 


v_L3 





epsilon.l + mu v.Ll [GeV"2] 
epsilon_2 + mu v_L2 [GeV"2] 
epsilon_3 + mu v_L3 [GeV"2] 



In case of the /iz/SSM the given are calculated from Equation (JSJ). The screen output 
moreover contains not only the one-loop on-shell neutralino and chargino masses, but also 
the tree level masses for comparison. Please note that for the correct calculation of the 
one-loop on-shell lepton masses a high accuracy of more than 20 digits is needed to ob- 
tain numerical stable results. In most cases for neutrino masses and the decay widths 15 
digits are sufficient, as long as all parameters are chosen real, implying that the diagonal- 
ization process isn't based on the squared mass matrices. To allow for a higher accuracy, 
adopt Selected_real_kind(15, 300) to (25,300) in /sphenooriginal/Control . F90 and com- 
pile CNNDecays. We note that a higher accuracy can result in significantly longer computing 
times. 
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